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Motivated by some previous works of Rudra et al we set to explore the background dynam- 
ics when dark energy in the form of New Variable Modified Chaplygin gas is coupled to dark 
matter with a suitable interaction in the universe described by brane cosmology. The main 
idea is to find out the efficiency of New variable modified Chaplygin gas to play the role of 
DE. As a result we resort to the technique of comparison with standard dark energy models. 
Here the RS11 brane model have been considered as the gravity theory. An interacting model 
is considered in order to search for a possible solution of the cosmic coincidence problem. A 
dynamical system analysis is performed because of the high complexity of the system . The 
statefinder parameters are also calculated to classify the dark energy model. Graphs and 
phase diagrams are drawn to study the variations of these parameters and get an insight into 
the effectiveness of the dark energy model. It is also seen that the background dynamics of 
New Variable Modified Chaplygin gas is consistent with the late cosmic acceleration. After 
performing an extensive mathematical analysis, we are able to constrain the parameters of 
new variable modified Chaplygin gas as m < n to produce the best possible results. Future 
singularities are studied and it is found that the model has a tendency to result in such 
singularities unlike the case of generalized cosmic Chaplygin gas. Our investigation leads us 
to the fact that New Variable Modified Chaplygin gas is not as effective as other Chaplygin 
gas models to play the role of dark energy. 
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I. INTRODUCTION 

Recent cosmic acceleration is a well-known and accepted fact in the cosmological society 
currently [l[ Qj. The root cause for this phenomenon is still under research. However of late there 
have been some speculations regarding the existence of a mysterious negative pressure component 
which violates the strong energy condition i.e. p + 3p < 0. Because of its invisible nature this 
energy component is aptly termed as dark energy (DE) 3j. 

Since the concept of DE flourished in the last decade, cosmologists all over the world started 
searching for a suitable model of DE. As a result various DE models have come into existence of 
late. DE represented by a scalar field 1 |4] is often called quintessence. Not only scalar field but 
also there are other Dark fluid models like Chaplygin gas which plays the role of DE very efficiently. 



As time passed extensive research was conducted and Chaplygin gas (CG) 



B,3, 



got modified into 

Generalized Chaplygin gas (GCG) [7Hlll] and then to Modified Chaplygin gas (MCG) In 
this context it is worth mentioning that dynamics of MCG in Braneworld was studied by Rudra 
et al [14f | . Other than these other forms of Chaplygin gas models have also been proposed such as 
Variable Modified Chaplygin gas (VMCG) [l5], New Variable Modified Chaplygin gas (NVMCG) 
[lfj], generalized cosmic Chaplygin gas fGCC G) fl?l| . Dynamics of GCCG in Loop Quantum 
cosmology was studied by Chowdhury et al in [18|]. The dynamics of GCCG in braneworld was 
studied by Rudra in Other existing forms of DE are phantom {^(J, k-essence 21], tachyonic 



field 



22], etc. 



The equation of state for NVMCG is given by, 

p = A(a)p-^, 0<a<l (1) 

Here we will consider A(a) = Aoa~ n and B{a) = Boa~ m , where Aq, Bq, a, m and n are positive 
constants. 

Currently, we live in a special epoch where the densities of DE and DM are comparable. The fact 
that they have evolved independently from different mass scales makes the fact more interesting. 
Given their non co-existence in evolution, comparable densities is quite an unexpected phenomenon. 
This is known as the famous cosmic coincidence problem. Till date several attempts have been 
made to find a solution to this problem j^ . A suitable interaction between DE and DM 

provides the best method of solution for this problem. It is obvious that a transition has occurred 



1 in the presence of a scalar field the transition from a universe filled with matter to an exponentially expanding 
universe is justified 
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from a matter dominated universe to dark energy dominated universe, by exchange of energy at 
an appropriate rate. Now in order to be consistent with the expansion history of the universe as 
confirmed by the supernovae and CMB data 25] the decay rate has to be fixed such that it is 



proportional to the present day Hubble parameter. Keeping the fact in mind cosmologists all over 
the world have studied and proposed a variety of interacting DE models 261431], 

Now, dark energy is not the only concept that can demonstrate the present day universe. The 
left hand side of the Einstein's field equation can also modified, to obtain suitable results. This 
modification however gives rise to the famous modified gravity theories, which in their own right 
can independently give us suitable models for our expanding universe. In this context Brane-gravity 
was introduced and brane cosmology was developed. A review on brane-gravity and its various 
applications with special attention to cosmology is available in 



32 



341 ] . In this work we consider a 



very popular model of brane gravity, namely the RS II brane. The main objective of this work is 
to examine the nature of the different physical parameters of the DE for the universe around the 
stable critical points in the brane model in presence of NVMCG. Effectiveness and success of the 
mathematical formulation of NVMCG will be studied. Impact of any future singularity caused by 
the DE in the brane world model will also be studied. 

This paper is organized as follows: Section 2 comprises of the analysis in RS II brane model. 
In section 3, a detailed graphical analysis for the phase plane is given. In section 4, some details 
regarding the mathematical construction of NVMCG is provided. In section 5 future singularities 
arising from the model are studied and finally the paper ends with some concluding remarks in 
section 6. 



II. MODEL 1: RS II BRANE MODEL 



Randall and Sundrum 
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361 ] proposed a bulk-brane model to explain the higher dimensional 
theory, popularly known as RS II brane model. According to this model we live in a four di- 
mensional world (called 3-brane, a domain wall) which is embedded in a 5D space time (bulk). 
All matter fields are confined in the brane whereas gravity can only propagate in the bulk. The 
consistency of this brane model with the expanding universe has given popularity to this model of 
late in the field of cosmology. 

In RS II model the effectiv e equat ions of motion on the 3-brane embedded in 5D bulk having 
^-symmetry are given by 



32, 



3j43g 



(4) G 



fJtV 



-A 4 g MI / + kIt^ + Kpl^ - E t 



(2) 



4 



where 



«I = g A«f , (3) 



A 4 = I 4(a 5 + i« ^2 . (1) 



2 ° V 6 
and 

and is the electric part of the 5D Weyl tensor. Here k 5 , A 5 , and A 4 are respectively the 
5D gravitational coupling constant, 5D cosmological constant, the brane tension (vacuum energy), 
brane energy- momentum tensor and effective 4D cosmological constant. The explicit form of the 
above modified Einstein equations in flat universe are 



3H" = A 4 + 4p+^ p 2 + -V ( 6 ) 
2 A Ak,a 



and 



2H + 3H* = A4-4P-^PP-^P 2 -^ (7) 

2A 2A A/c| 

The dark radiation IA obeys 

U + AHU = (8) 

where p = p n vmcg + Pm and p = p n vmcg + Pm are the total energy density and pressure respectively. 

As in the present problem the interaction between DE and pressureless DM has been taken into 
account for interacting DE and DM the energy balance equation will be 

Pnvmcg + 3.ff (1 + uj nvmcg ) p nvmC g = ~Q, for NVM CG and (9) 

p m + 3Hp m = Q, for the DM interacting with NVMCG. (10) 
where Q = 2>bHp is the interaction term, b is the coupling parameter (or transfer strength) and 

P — Pnvmcg Pm IS t 



p + 3H (p + p) = 



le total cosmic energy density which satisfies the energy conservation equation 
IB. 
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Since we lack information about the fact, how does DE and DM interact so we are not able to 
estimate the interaction term from the first principles. However, the negativity of Q immediately 
implies the possibility of having negative DE in the early universe which is overruled by to the 
necessity of the second law of thermodynamics to be held 4l[. Hence Q must be positive and 



small. From the observational data of 182 Gold type la supernova samples, CMB data from the 
three year WMAP survey and the baryonic acoustic oscillations from the Sloan Digital Sky Survey, 
it is estimated that the coupling parameter between DM and DE must be a small positive value 
(of the order of unity) , which satisfies the requirement for solving the cosmic coincidence problem 
and the second law of thermodynamics j^ . Due to the underlying interaction, the beginning of 
the accelerated expansion is shifted to higher redshifts. The continuity equations for dark energy 
and dark matter are given in equations (9) and (10). Now we shall study the dynamical system 
assuming A4 = U = (in absence of cosmological constant and dark radiation). 

A. DYNAMICAL SYSTEM ANALYSIS 

In this subsection we plan to analyze the dynamical system. For that firstly we convert the 
physical parameters into some dimensionless form, given by 

, Pnvmcq Pm Q> 

* = lna, - = ^T, - = ^ V = ^p (11) 

where the present value of the scale factor ao = 1 is assumed. Using eqns. (1), (6), (7), (9) and 
(10) into (fTT|) we get the parameter gradients as 



g = -36 (u + v) - 3 (u + uuifjjj) - 6§u (12) 



^ = 3b(u + v) - 3v - Q^pv (13) 



and 



d I=y{ i -4) «") 

where uin^mcg is the EoS parameter for NVMCG determined as 

(RSII) = Pnvmcg 1 (._ By n ~ m \ 

HVmC9 Pnvmcg y n X^ \ u a+l X m-n+a + l ) > ^ 
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H 



u + v 



U + V 



and 



nvmcg ^ «2 („ + „) 



-1 (u( W ^ + 2)+2«) 



2A 

u + v 



K 2 (u + V) 



(16) 



(17) 



1. CRITICAL POINTS 

The critical points of the above system are obtained by putting ^ = ^ = ^ = 0. But due to 
the complexity of these equations, it is not possible to find a solution in terms of all the involved 
parameters. So we find a solution for the above system, by putting the following numerical values 
to some of the parameters appearing in the system. We take, 

a = 0.5, b = 0.5, n = 1, m = 1 

and obtain the following critical point, 



0.3125 



0.1875 



1.02337 x 10~ 20 / 8.19379 x 10 18 B\% 1.14512 x 10 19 ylA 4 



U r . = 



V r = 



K 



K 



A 5 



K 



K 



(18) 



The critical point correspond to the era dominated by DM and NVMCG type DE. For the 
critical point (u c ,v c ), the equation of state parameter given by equation (15) of the interacting DE 
takes the form 



Co' 



(RSII) _ Pnvmcg 



1 



nvmcg 



Pnvmcg Vc-^- 71 



A 



By? 



u a+l xm - n +a+l I ' 



where 



X = 



2A 



u c + v c 



K 2 (u c + V c ) 



(19) 



(20) 



Fig 1 : The dimensionless density parameters u, v and y are plotted against each other in a 3D- 
scenario. Other parameters are fixed at a = 0.5, b = 0.5, Ao = 1/3, Bq = 3, n = 2, m = 1, A = 1.5, 
and k = 0.2. 



Fig 2 : The dimensionless density parameters are plotted against e-folding time. The initial 
conditions are v(lA) = 0.05,u(l.l) = 2.5 and y(l.l) = 2.8. Other parameters are fixed at 
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Fig.2 Fig.3 Fig.4 



a = 0.5, b = 0.001, A = 1/3, B = 3, n = 5, m = 1, A = 1.5, and k = 0.2. 

Fig 3 : The dimensionless density parameters are plotted against e-folding time. The initial 
conditions are v(l.l) = 0.05,n(l.l) = 2.5 and y(l.l) = 2.8. Other parameters are fixed at 
a = 0.5,6 = 0.1, A) = 1/3, S = 3, n = 5, m = 1, A = 1.5, and k = 0.2. 

Fig 4 : The dimensionless density parameters are plotted against e-folding time. The ini- 
tial conditions are v(l.l) = 0.05,n(l.l) = 2.5 and y(l.l) = 2.8. Other parameters are fixed at 
a = 0.5, b = 0.5, Aq = 1/3, Bq = 3, n = 5, m = 1, A = 1.5, and k = 0.2. 
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Fig 5 : The phase diagram of the parameters u(t) and v(t) depicting an attrac- 
ted solution. The initial conditions chosen are v(l) = 0.05, u(0) = 2.5, 7/(1) = 1.8 
(green); u(l) = 0.06,7/(1) = 2.6,7/(1) = 1.9 (blue); u(l) = 0.07,7/(1) = 2.7, y(l) = 2 
(red); v(l) = 0.08, u(l) = 2.8,7/(1) = 2.1 (gold). Other parameters are fixed at 
a = 0.5, b = 0.01, A = 1/3, B = 3, n = 5, m = 1, A = 1.5 and k = 0.2. 

Fig 6 : The dimensionless density parameters u(t) and y(t) are plotted against e-folding 
time. The initial conditions are v(l.l) = 0.05,^(1.1) = 2.5 and 7/(1.1) = 2.8. Other parameters 
are fixed at a = 0.5, b = 0.01, Aq = 1/3, Bq = 3, n = 5, m = 1, A = 1.5, and k = 0.2. 
Fig 7 : The dimensionless density parameters v(t) and y{t) are plotted against e-folding time. The 
initial conditions are u(l.l) = 0.05,n(l.l) = 2.5 and 7/(1-1) = 2.8. Other parameters are fixed at 
a = 0.5, 6 = 0.001, A = 1/3, B = 3, n = 5, m = 1, A = 1.5, and « = 0.2. 

Fig 8 : A 3D-phase portrait of the dimensionless density parameters u(t), v(t) and y(t) is 
plotted against each other. The initial conditions are v(l.l) = 0.05,7x(l.l) = 2.5 and 7/(1.1) = 2.8. 
Other parameters are fixed at a = 0.5,6 = 0.001, Aq = 1/3, Bo = 3, n = 5, m = 1, A = 1.5, and 
k = 0.2. 
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Fig.9 

Fig 9 : The deceleration parameter is plotted against the EoS parameter. Other parameters 
are fixed at a = 0.5, b = 0.001, Aq = 1/3, Bq = 3, n = 5, m = 1, A = 1.5, and k = 0.2. 



Fig 10 : The statefinder parameter r is plotted against the EoS parameter. Other parameters 
are fixed at a = 0.5, b = 0.001, Aq = 1/3, Bq = 3,n = 5, m = 1, A = 1.5, and k = 0.2. 
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Fig 11 : The statefinder parameter s is plotted against the EoS parameter. Other parameters are 
fixed at a = 0.5,6 = 0.001, A) = 1/3, B = 3, n = 5, m = 1, A = 1.5, and k = 0.2. 

Fig 12 : The ratio of density parameters is shown against e-folding time. The ini- 
tial conditions chosen are v(l)=0.05, u(l)=2.5, y(l)=1.8. Other parameters are fixed at 
a = 0.5, b = 0.001, A = 1/3, B = 3, n = 5,m = 1, A = 1.5, and k = 0.2. 
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Fig. 12 



2. STABILITY AROUND CRITICAL POINT 

Now we check the stability of the dynamical system (eqs. (12) and (13) and (14)) about the 
critical point. In order to do this, we linearize the governing equations about the critical point i.e., 

u = u c + Su, v = v c + 5v, y = y c + 5y (21) 

Now if we assume / = ^ , g = and h = then we may obtain 

5 (£) = [duf]c 6u + [dvf]c 5v + [dyf]c 5y (22) 



S ( ^\ = [d u g] c 5u + [d v g] c 5v + [d y g] c 5y (23) 



and 



5 (S) = [duk]JU + [dvk]c 5V + [dyk]c 6V (24) 



where 
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d u f = -3 



2-l-m-n-a u -l-a y -m-n / \(- U - V + 4j) 



A(« + v) 2 (-l + (u + v)k 2 ) 2 \ (u + v) 2 



-2 n Bv(u + v) 2 y n K 2 (u + 2mu + ua - va - (u + v){mu - va) K 2 ) " ^r^ l + 2 1+m+a u a y m 



/ A(— U — V + \ i , 

x(-1 + (u + v)k 2 )X — 1 ' " • 1 1 ' ' - 



(Aw(-2nu - v + (u + v)(nu + v)k 2 ) + 2 n (u + u)V 



X (_l + ( U + „)« 2 )(-l + 6 + 2(2^ + t,)/, 2 ) ( A( ( 1 + ^2^ ) 



-m—n—a\ 



(25) 



^/ = -3 



2-l-m-n-a u -a y -m-n i \(- U - v + ' 



-m— n— « 



A(lt + v) 2 (-l + (it + 1>)k 2 ) 2 I (li + 77) 2 



X 



no o I X( — U — V + 4t) 

2 n 5(u + v) 2 y n K 2 {-u - 2v{l + m + a) + (it + u)(u + v(l + m + a))rc 2 ) I — - K 



+ 2 1+m+a u a 



(Av(-it - 2m; + (u + v)(u + m;)/c 2 ) + 2 n {u + w) z 



x(6 + 2^ 2 )(-l + (u + t>) K 2 ) f A( " 

\ (u + v) z 



(26) 



2- n 3Anuvy 



dyf = 



-i-n f A( M ' 2- 1 - m - a 35mr tt «( ?J + v)y- 1 - m K 2 ( H ( " + ^ } 



+ 



(-1 + (u + u)k 2 )A 



(27) 



<9„fif = 3 



2 -l-m-n-a u -l-a y -m-n / \(-u - V + 4j) % 



-m— ra— « 



A(u + v) 2 {-\ + (u + v)k 2 ) 2 I (u + vf 



—2 n Bv(u + v) y n K (u + 2mu + ua — va — (u + v){mu — va 



/A(-u-,+^y +2l+m+v+ ^ 
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/ l \ \ m ~\~ a 



(it + v) s 



(Au(-2mt + + + i;)k 2 ) + 2 n (u + u)V 



x(6-2i;k 2 )(-1 + (u + wK) 



2 / A(-«-u + 4r) 



(it + w) 5 



(28) 



S„5 = 3 



A(« + v) 2 {-\ + {u + v)k 2 ) 2 \ (u + v) 2 



2 n B(u + v) 2 y n K 2 {-u - 2v{\ + m + a) + (it + v)(u + v(l + m + a))/t 2 ) f — " + ^ ^ 



2l+m+a u a 



( \ I i s \ m+a 



(u + w) 2 



(Au(u + 2m; - (u + v)(u + uv)k 2 ) + 2 n (u + v) 2 y r ' 



x (-1 - b + 2(u + 2?;)k 2 )(-1 + (u + f )k 2 



'A(-u-u + ^) x 



u + v) 



(29) 



2-" :3,1»» rir 1 ' 11 (^^M) 2- 1 - m - a Bmu- a v(u + v)y- 1 ~ m K 2 (^^1 



d y g 



i . \ — m— a 



tt + f 



(-1 + (u + v)k 2 )A 



(30) 



a,/i = -3 



2 -l-m-n-a u -l-a y -m-n I \(- u - v + J^) \ 



A(« + w) 2 (-l + (u + v)k 2 ) 2 I (u + v) 2 



X 



-2 n £(u + v) 2 y n K 2 (-u - 2mu - ua + va + (« + u)(mu - to)k 2 ) [ ^— ^ — - % — | + 2 1+m+c V +a y m 



x(-1 + (u + w)k 2 )A 



~A(-n-t; + ^)' 
(u + w) 2 



m+a 



(-A(2nu + v)+ A(u + «)(nu + v)k 2 - 2 1+n (u + v) 2 y n n 2 



x(— 1 + (u + 



(u + f)" 



(31) 
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d v h = —3 



2 -l~rn-n-a u -a y l-m-n / \(- U - V + 4j) ' 



-m— n— « 



A(u + v) 2 (-l + {u + v)k 2 ) 2 \ (u + v) 2 



-2 n B(u + v) 2 y n K 2 {-l -2m-2a + {u + v)(m + o)k 2 ) ( ^ U V + 2 1+m+a u a y m 



(U + V) 



x(-l + (u + v)k 2 )X 



(u + v) 2 



m+a 



(Au(-1 + 2n - (n - l)(u + v)k 2 ) + 2 l+n {u + i>) 2 yV 



x(— 1 + (u + 



(u + V ) 2 



(32) 



<9„/i 



(u + u)(-l + (u + t;) k 2 )A 



, / (-u - v + -4) A 

2 — 1 — m—n—oiy—oiy—m— n | \ 



(it + f) z 



-3 x 2™5 (m - 1) (u + w) 2 yV ^ " " + f^ A > ) - 2 1+m+Q ?/*y m (-1 + (u + u) k 2 ) Ax 

V (u + v) I 



{ - U - V + ^ X ) m+a {3A(n-l)u + 2^ U + v)y-(6(u + v) K 2 -7) ( { ~ U ~ V + & V 



(it + w) 2 y i 

The Jacobian matrix of the above system is given by, 



(u + vY 



J, 



(RSI I) 
(u,v) 



( 81 

Su 


§1 

Sv 


§1 \ 

Sy 


Sg 
Su 


Sg 
Sv 


Sjl 
Sy 


Sh 
\ Su 


Sh 
Sv 


Sh 

Sy I 



The eigen values of the above matrix are calculated at the critical point (u c ,v c ) and are found 
to be Ai = 5.99118, A 2 = -3, A 3 = -2.67056. Hence it is a Saddle point. 



3. NATURE OF COSMOLOGICAL PARAMETERS 



Deceleration Parameter: 
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In this RSII model, the deceleration parameter q can be obtained as 

\ (, (RSII) 
(RSII) _ , 3 1 2A [Vnvmcg — I) ~ \V + W nvm cg — j 

" ~ ^ (1+^) (33) 

which can be written in terms of dimensionless density parameter 0, nvmcg = Pnv ™ C3 a s in the 
following 

3 { ("2A + "0 w nvmcJ ^nvmcg + (l + f ) } 

9 < »"> = -l + 5 i ^ i (34) 

Now since fl nvmcg = Pnv ™ C3 = ^ and assuming § = e (i j 5//) we get, 

.,*,,„ .^ {(i-^^'^+a+^n,); (35) 

Considering only the first stable critical point, such that (u,v) — > {ui c ,v\ c ), using (f35l) we get, 

(36) 

If ^(.RS//) = — jz — (rsii)\ ^ — > x (RSii) = 0, we have q = -1, which confirms the 



{(l_Ssp> 


(i?577) tt /- n 
Wnvmcg ^ + I 1 + e (RSII) / 


} 







T7T (RSii)\ ~7 

i Z IMnvmcg J ^ 1 c ~r 1 c 

accelerated expansion of the universe. When emsil) = ~ 2 we have q = —oo. Therefore we have 
super accelerated expansion of the universe. 



In this scenario, the Hubble parameter can be obtained as, 



H = vT—t (37) 

^(RSIIY 



where the integration constant has been ignored. Integration of (|37p yields 



(t) = a t 3X (Rsn) (38) 

which gives the power law form of expansion of the universe. In order to have an accelerated 
expansion of universe in RSII brane we must have < -X^st/) < §• Using this range of X^ R gjj^ 
in the equation q^ sn = — 1 + §-X7RSjn, we get the range of qi RSH ^ as — 1 < q c RSH ^ < 0. This is 
again consistent with the accelerated expansion of the universe. 
2. Statefinder Parameters 

As so many cosmolo gies ! models have been developed, so for discrimination between these 



contenders, Sahni et al 



431 ] proposed a new geometrical diagnostic named the statefinder pair 
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{r,s}. The statefinder parameters are defined as follows, 

d r — 1 



(39) 



aH 3 ' 3(g-l/2) 
where a is the scale factor of the universe, H is the Hubble parameter, a dot denotes differentiation 
with respect to the cosmic time t. 

Here we calculate the statefinder parameters {r, s} in order to get relevant information of DE 
and DM in the context of background geometry only without depending on the theory of gravity. 
The expressions of the statefinder pair eq. ()39p in the RSII model can be obtained in the form 

r( R sn) = (l " (1 - 3X {RSII) ) . (40) 

and 

S (RSII) = X (RSII) (41) 

III. GRAPHICAL ANALYSIS 

Graphs are obtained and phase diagrams are drawn in order to determine the type of critical 
point obtained in this model. Below we discuss the results obtained in detail: 

The dimensionless density parameters u and v and y are plotted against each other in figure 1. 
From the figure we see that v decreases, and u increases during evolution of the universe. In figs. 
2,3 and 4, the density parameters are plotted against time. It is evident from the figures that the 
density of DM decreases while the density of DE increases as the universe evolves with an increase 
in scale factor a. So this result is consistent with the well known idea of an energy dominated 
universe. In these figures it is seen that with an increase in the value of interaction the values 
of u and v become more and more comparable to each other, which is quite an expected result. 
So we have a possible solution of the cosmic coincidence problem. A comparative study between 



14, [3, Q and the current, work reveals thai DE domination over DM i> much less pronounced 
in case of NVMCG than in case of other Chaplygin gas models like, MCG or GCCG. This is a 
very interesting feature in the character of NVMCG. Moreover a comparative study of the current 
work with [l8] reveals the fact that density of DM is much more comparable to the density of 
GCCG than to that of NVMCG. Hence it can be speculated that NVMCG is perhaps less 
effective to play the role of DE in comparison to MCG and GCCG as well. 

Figs. 5, 6 and 7 shows the phase portrait of the density parameters of DE and DM. In fig. 
5, we see a phase diagram between density parameters u(t) and v(t). Figs. 6 and 7 shows the 
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phase diagram between u(t), y(t) and y(t), v(t) respectively. Finally in fig. 8, a 3D-phase diagram 
between all the three density parameters is obtained in a single system. As already stated before 
that the critical point obtained in this system is a Saddle point and hence there always remains 
a question on the stability of the system. In fig.9, a plot of deceleration parameter, q is obtained 
against the EoS parameter, oj. It is seen that q remains in the negative level thus confirming the 
recent cosmic acceleration. Figs. 10 and 11 show the plot of the statefinder parameters r and s 
respectively against the EoS parameter, oj. It is known that in case of ACDM model r = 1 and 
s = 0. Here we see that in fig. 11 except at oo = 0, the values of s correspond to that of the ACDM 
model. In fig. 10, we see that the values of r is quite different from 1 corresponding to the values 
of oo when s = in fig. 11. 

This gives the deviation of the model from the ACDM model. Finally in fig. 12, the ratio v/uis 
plotted against x = In a. The decreasing trajectory confirms the existence of an energy dominated 
universe with progressive values of scale factor. 

IV. SOME NOTES ON THE MATHEMATICAL CONSTRUCTION OF NEW 
VARIABLE MODIFIED CHAPLYGIN GAS 

We know that NVMCG is basically an extension of modified Chaplygin gas, and its mathe- 
matical formulation is based on the modification of the EoS of MCG in a way that makes it more 
suitable as a candidate to play the role of DE. Our motive is to investigate how far successful is 
NVMCG over MCG as a DE and consequently how far essential or justified is this modification. 
So in this section we consider modified Chaplygin gas as a standard model of dark energy and 
perform a comparative study with NVMCG. The EoS for MCG is given by, 



P — Ap mcg — 

rmcg 

where A, B and a are positive constants. For negative pressure we should have 



(42) 




(43) 



Now to get negative pressure in case of NVMCG we should have, 




(44) 



Case I: for m = n 



From relation (43), we have, p nvrnC g < 
MCG for B = B and A = A. 




a + l 



. The above value of p nvrnC g coincides with that of 
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Case II: for m> n 

Let m — n = mi, where mi > 0. Therefore relation (43) becomes, p < (%) a=+r. Now since 

a > 1 in an accelerating universe, so a«+i > 1. Let a^+i = 7712. So we get p n vmcg < ( a^J ° m 2> 
Now since 777,2 > 1, p m cg < Pnvmcg if Aq = A and -Bo = B. Here it can be seen that the density of 
NVMCG increases due to the introduction of scale factor, o in its EoS. It can be seen from the EoS 
of NVMCG that this effect increases the magnitude of the first term and decreases the magnitude 
of the second term, which eventually increases the the value of pressure as a whole. This push 
towards the towards the positive region in the value of pressure can be speculated as a basic flaw 
in the mathematical construction of NVMCG, since it reduces its efficiency as a DE. Hence m > n 
is not at all acceptable as far as the concept of dark energy is concerned. 
Case III: for m < n 

In this case m — n = —mi. Using this we get, p < - _ mi — . For an accelerating universe 

a "+ 1 

a > 1, which implies a^+i > 1. Therefore p < ^ ° +1 • Now since m<i > 1, which implies 
Pmcg > Pnvmcg, if ^-o = ^4 and Bq = B. So there is a relative decrease in the value of density of 
NVMCG compared to its counterpart MCG. This in fact pushes the value of pressure towards the 
more negative region, thus enhancing the ability of NVMCG to play the role of DE. This is the 
most acceptable case and the EoS of NVMCG should be constrained by m < n. 



V. STUDY OF FUTURE SINGULARITIES 



We know that any energy dominated model of the universe will result in a future singularity. 
As a result the study of dynamical model of universe in the presence of DE and DM is in fact 
incomplete without the study of these singularities, which are the ultimate fate of the universe. It 
is known that the universe dominated by phantom energy ends with a future singularity known as 



Big Rip 44], due to the violation of dominant energy condition (DEC). But other than this there 



are other types of singularities as well. Nojiri et al [4| studied the various types of singularities 
that can result from an phantom energy dominated universe. These possible singularities are 
characterized by the growth of energy and curvature at the time of occurrence of the singularity. 
It is found that near the singularity quantum effects becomes very dominant which may alleviate 
or even prevent these singularities. So it is extremely necessary to study these singularities and 
classify them accordingly so that we can search for methods to eliminate them. The appearance 
of all four types of future singularities in coupled fluid dark energy, F(R) theory, modified Gauss- 
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Bonnet gravity and modified F(R) Horava-Lifshitz gravity was demonstrated in [20j]. The universal 
procedure for resolving such singularities that may lead to bad phenomenological consequences was 
proposed. In Rudra et al JjJ] it has been shown that in case of Modified Chaplygin gas(MCG), 
both Type I and Type II singularities are possible. However in [l9] it was shown that GCCG does 
not result in any type of future singularity. 



A. TYPE I Singularity (Big Rip singularity) 

If p — > oo , \p\ — > oo when a — > oo and t — > t s . Then the singularity formed is said to be the 
Type I singularity. 

In the present case by considering the NVMCG equation of state from equation (1) we see that 
when a — > oo, \p\ — > 0. Therefore we see that there is no possibility for Typel singularity, i.e., Big 
Rip singularity, in case of NVMCG. 



B. TYPE II Singularity (Sudden singularity) 

If p — > p s and p s ~ 0, then \p\ — > — oo for t — > t s and a — > a s , then the resulting singularity is 
called the Type II singularity. 

Considering the equation of state for NVMCG, We see that if p s ~ 0, then \p\ — > — oo for t — > t s 
and a — > a s . Hence there is a strong possibility of the type II singularity or the sudden singularity 
in case of NVMCG. 



C. TYPE III Singularity 

For t — > t s , a ^ a s , p — >oo and \p\ — > oo. Then the resulting singularity is Type III singularity. 
It is quite evident from the equation of state of NVMCG that it supports this type of singularity. 



D. TYPE IV Singularity 

For t — >• t s , a — > a s , p — >• and \p\ — >• 0. Then the resulting singularity is Type IV singularity. 
Investigation shows that this type of singularity is not supported by NVMCG type DE. 

As a remark, one should stress that our consideration is totally classical. Nevertheless, it is 
expected that quantum gravity effects may play significant role near the singularity. It is clear 
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that such effects may contribute to the singularity occurrence or removal too. Unfortunately, due 
to the absence of a complete quantum gravity theory only preliminary estimations may be done. 

VI. CONCLUSION 

In this work we have considered New variable modified Chaplygin gas and tried to determine its 
efficiency to play the role of dark energy in an universe described by RSII brane. A numerical system 
study was carried out in order to throw some light on the dynamics of the dark energy. The system 
was formed and a solution was obtained. An eigen value analysis of the system at the critical point 
showed that the system was far from attaining stability since it produced a saddle point. Plots were 
obtained to get a clear idea about the result of our analysis in a both qualitative and quantitative 
aspect. It was found that NNMCG in RSII brane model is perfectly consistent with the idea of an 
energy dominated universe, but the DE domination over DM is much less pronounced in case of 
NVMCG when compared to that with Modified Chaplygin gas or generalized cosmic Chaplygin gas. 
This is an important result indeed. It was also discovered from the plots that with the increase in 
the magnitude of interaction the values of DE and DM became more and more comparable to each 
other thus providing a solution of the cosmic coincidence problem at higher interaction scales. The 
plot of the deceleration parameter revealed that the model is perfectly consistent with the notion 
of an energy dominated universe. The unique trajectories in the plots of statefinder parameters 
differentiated the model from other DE models. An extensive study regarding the mathematical 
formulation of NVMCG was performed, and it was found that m < n best suited the nature of 
NVMCG as a DE. Hence we have been able to constrain the parameters of NVMCG 
to give the best possible results. Finally the future singularities were studied and the model 
was found to be affected by some of them, quite unlike the generalized cosmic Chaplygin gas. In a 
nutshell it can be said that the performance of NVMCG as a dark energy is quite moderate. When 
compared with other DE models like MCG or GCCG, the weakness of NVMCG is quite visible. 
Hence as a conclusion we speculate that NVMCG is not quite the brightest contender to play the 
role of DE and there is a lot of scope for improvement as far as the mathematical aspect of the 
model is concerned. 
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